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Interaction Protected Topological Insulators with Time Reversal Symmetry 
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Anderson’s localization on the edge of two dimensional time reversal (TR) topological insulator 
(TI) is studied. For the non-interacting case the topological protection acts accordingly to the Z2 
classification, leading to conducting and insulating phases for odd and even fillings respectively. 

In the presence of repulsive interaction the phase diagram is notably changed. We show that for 
sufficiently strong values of the interaction the zero temperature fixed point of the TI is conducting, 
including the case of even fillings. We compute the boundaries of the conducting phase for various 
Hllings and types of disorder. 


I. INTRODUCTION 


Time reversal non-interacting TIs are realized in ma¬ 
terials with strong spin orbit interaction in twcP^^ and 
thre^^lt^ dimensions. Also known as quantum spin Hall 
insulators, these materials have an insulating bulk, while 
hosting gapless surface states. 

For the non-interactingcase the classification of dis¬ 
ordered TIs is completeSHHl. In two dimensions, non in¬ 
teracting TR invariant TIs are classified by a Z 2 topo- 
logic al inva riant, according to the number of helical edge 
stateJi^E^. The back scattering by TR invariant disor¬ 
der is possible only between states that are not Kramer’s 
partners. Therefore for an odd number of the helical 
states the conducting state is protected. For an even 
number of helical edge modes the electrons can be local¬ 
ized completely by scattering among non-Kramers pairs. 
This state is therefore equivalent to a trivial insulator in 
agreement with a more general Haldane criteriorP^. 

Interacting integer and fractional TI’s are a subject of 
active research. Interaction may lead to strongly corre¬ 
lated ground state^i^lt2ll with fractional excitations and 
non trivial statistic^22H24]_ classification of TI’s in the 
interacting case is yet unknown. An approach, based on 
thermal response and its relation to a quantum anomaly, 
valid beyond single particle picture, was proposed in 
Refill. 

It is commonly accepted, that the charge transport in 
the ideal TI occurs via protected a single helical edge 
mode, with the universal quantized conductance 2e^//i. 
In reality the conductance differ from this value due to 
back scattering processes. The latter may occur via the 
combination o f the two electron scattering and the dis¬ 
order potentiaPSIHI^ coupling to the bulk via electron 
puddled or due to the magnetic impurities. In the later 
case the interaction stabilizes the conducting phase, and 
the quantum phase transition as function of Luttinger liq¬ 
uid (LL) parameter is predicted for the fractional TlP^I 

In this work we study the localization by TR disorder 
on the edge states of a TI in the presence of repulsive 
interaction. Although we focus on the TI’s at integer 
fillings {ly), that in the absence of the disorder possess ly 
helical edge states, the same analysis applies for narrow 
stripe of TI at = I.^ The inclusion of TR disorder 


drives the non-interacting system to a state with a single 
or no helical edge states. We show that the presence 
of the repulsive interaction can stabilize the conducting 
phase. 

We model the disorder by a short range static poten¬ 
tial that due to the spin orbit interaction mixes different 
helical states, except those that are connected by the TR 
symmetry. We consider a generic finite range interaction 
between the electrons, with all possible matrix elements 
allowed by symmetry. We consider the case of a sin¬ 
gle impurity and the random disorder, with a scattering 
length shorter that the sample size. We perform one loop 
renormalization group (RG) anaWsis, analogous to Kane- 
Fisheil^and Giamarchi-Schultd^ study of localization in 
one dimensional systems. 

Our analysis shows that the low energy fixed point 
is determined by the magnitude of the interaction and 
its effective radius. For interaction stronger that some 
critical value the low temperature phase is conducting. 


II. THE MODEL 


The appearance of the helical edge states can be un¬ 
derstood on the level of non interacting electrons. In 
the presence of a Rashba spin-orbit (SO) interaction the 
single particle Hamiltonian is given by 

H = ^ oisoip X a) ■ \/V{x,y) + V{x,y), (1) 

where rrie is effective mass of an electron and aso is the 
strength of the SO coupling. For the parabolic potential 
V{y) = y'^ j2meagQ, the Hamiltonian Q 
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( 2 ) 


corresponds to two replicas of fermions subject to oppo¬ 
site magnetic fields \B\ = a^Q. For the integer fillings, 
iz = a'^QA/i^Q (A being an area of sample and <i)o = hcje 
the flux quantum) the bulk forms an incompressible state 
with a gap of size ha'^Q/rrie- 

An addition of a smooth confining potential curves the 
Landau levels, as shown in Figg leading to v gapless 
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helical edge state^i^*^^. Assuming that the single par¬ 
ticle gap formed in the bulk is not closed the effects of 
interaction can be taken into account within the helical 
edge states. This phenomenological approach can be mi- 
crosc opical ly justified within the sliding Luttinger liquid 
modepSEB, However, the resulting helical edge descrip¬ 
tion is believed to be the correct low energy model, valid 
beyond the sliding LL approximation. 

To account for the interaction, it is natural to 
pass to the bosonic description, defining bosonic 
fields Ri/Li, related to the right/left density compo¬ 
nents by pR^i = dxRi/2'K and pL,i = -dxLi/2^. 
These fields satisfy the canonic commutation relations 
[Rj{x),Rj>{y)] = -[L,{x),Lj,{y)\ = iTTS^{x - y)5jp. 
The electronic operators are represented as i^Rj = 
j^2'Ka, and j'J2 tui, with a a short 

distance cut-off. 

In the absence of the Umklapp and 2kp electron- 
electron scattering the interaction between v modes, con¬ 
sistent with TR symmetry, is represented by the following 
action 

S' = ^ y dxdt (d^^^JCdt^ - , (3) 

where we use the compact notations $ = {R,L), R = 
(i?i, i? 2 , ■ • •, Sj/) and similarly for L. Here the matrix 
JC encodes the commutation relations of the fields and 
can be written as /C = CTz ® Ii,, where the tJz is a Pauli 
matrix that acts in the right/left movers subspace while 
Ii/ is the identity matrix in the space of v modes. The 
positive definite matrix Ai accounts for interaction. The 
helical edge modes are separated in space by a distance 
d. 

The symmetry under TR requires that {T,/C} = 
[T, A4] = 0, where T = CTj, ® Ii, is a time reversal symme¬ 
try operator. This restricts the interaction matrix A4 to 
the form 


M = 


Mfw Mbw 
Mbw Alfw 


I 2 0 MfwMbw, (4) 


where describes the forward interaction between 

the copropagating modes pR^i and pRj (similarly for left 
movers). is a an element of a symmetric matrix 

in the channel space that describes the backward interac¬ 
tion between Ri and Lj. We assume that the interaction 
between helical modes i and j is translationally invariant, 
and depends only on the relative distance \i — j\. In the 
absence of disorder the spectrum of this model is gapless. 

The presence of impurities may dramatically change 
the states of edge modes. TR invariant disorder mixes 
helical states that belong to different Kramer’s pairs and 
induces backward scattering processes. We consider two 
cases: (a) single impurity scattering and (b) random dis¬ 
order. Single impurity scattering is the dominant process 
if the mean free path of the electrons is larger than the 
sample size. In the opposite case (the electrons’ mean 
free path is smaller than the sample size) the localization 
is dominated by multiple scattering. 



FIG. 1. (color on line) Helical edge modes. The wavy lines 
represent the back-scattering events between non-Kramers 
pairs (in red and blue). Insert.- Schematic band structure 
of a TI. 


A. Single Impurity 

We analyze the single impurity case first. The sin¬ 
gle impurity, located on the edge backscatter between 
states i and j that are not connected by TR symmetry, 
~ i( 0 )V^L,i( 0 ), where pij is proportional to 

the impurity potential ai k = The renormalization 
of the strength pij is a straightforward generalization of 
standard Kane-Fisher analysi^^ 

^ = (1-A,,)m,,. (5) 

Here A^- is the scaling dimension of the scattering pro¬ 
cess (V’L(^)V'L.i('r)^Ij(0)r/'fl,*(0)) - From 

Eq. (§ it follows that a single impurity is an irrelevant 
perturbation if Ay > 1 . 

The scaling dimension A is controlled by the back¬ 
ward interaction matrix Mbw. For the simple case v = 2 
only two helical modes propagate on the edge. If the 
separation between the helical states is larger than the 
interaction radius the effective Hamiltonian is given by 
R = i/fw + Mbw with 

Mfw = ^ mR^f + 

i—1 ^ / 

0 ^ r 

Mbw = “% X! / dxdxRidxLi. (6) 

^ i=i •' 

Here vr the Fermi velocity while 54 and param¬ 
eterize the forward and backward interactions respec¬ 
tively. Defining the fields pi = {Ri — Li)f^/2 and 
9i = {Ri+Li)/'/2, the above Hamiltonian can be written 
in a Luttinger liquid (LL) form 

. ( 7 ) 
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where u = (z;f + 54/'^)\/l ~ (-^2)^ is the renormalized 
sound velocity, A 2 = and K = is the 

LL parameter [K < 1 for repulsive and K > 1 for at¬ 
tractive interaction within the mode). In this case the 
single impurity operator has the scaling dimension 
Ai 2 = i^/2 -I- 1/2K^ so the single impurity is irrelevant 
for any interaction. This is in stark contrast with the 
non helical LL where the scaling dimension of the disor¬ 
der operator is K, so the disorder is irrelevant only for 
an attractive interactiorP^. If the LL’s are different the 
scaling dimension is controlled by both LL parameters 
Ai 2 = {Ki + 1/Ki + K 2 -\- 1 /A' 2 )/ 4 . This result implies 
that the scattering between different edge states is ac¬ 
companied by the zero bias anomaly that suppresses the 
probability of this process. This is in contrast to Kane- 
Fisher cas^^, where the back scattering occurs in the 
same LL and has no zero bias anomaly suppression. For 
TI this processes is forbidden by the TR symmetry. 

For the case where the helical edge states are located 
within the radius of interaction (or for any long range 
interaction potential) the interaction matrices Q are 




fvF+gl gl \ 

V gl VF + gl) 


ALbw 


5 ° gl 
gl gl 


( 8 ) 


The scaling dimension A 12 is 


We now follow the steps of Giamarchi-Schultz renor¬ 
malization group analysid^. For the weak disorder one 

findPI 


dW,, 

d£ 


(3 - 2A,,)W„-, 


(14) 


where A ^,- = A|j_j| is the scaling dimension of scattering 
process ( 12 ) between helical states i and j allowed by 
TR symmetry. In the conducting phase the disorder is 
an irrelevant perturbation, and all Wij flow to zero. This 
requires A|j_^| > | for all pairs i,j. Let us consider two 
limiting cases: (i) disorder that mixes only the nearest 
modes Wij ^ ITdij+i; (ii) the disorder that mixes the 
modes uniformly Wij ^ W. All physical realizations lie 
in between these two limits. 


The simplest situation is realized for v = 2 where 
the limits (i) and (ii) coincide. In that case, the scal¬ 
ing of the disorder operator is given by Eq.(|9). In the 
absence of inter-mode interaction the scaling dimension 
of a back scattering operator is A 12 = {K + K-^)/2. 
Therefore the system flows to a conducting fixed point 
for AT < (3 — •\/5)/2. In the presence of inter mode inter- 


Ai 2 = - {Fi 




where, for repulsive interactions 


Fi 

F 2 


1 + -^4 ~ ^2 ~ ^^2 

i + xl + xl + xV 


'l-Xl + X°-Xl 
1 — Al — Ao -I- An ’ 


Fi G [0,1] 


F 2 G [0, 00 ] 


( 9 ) 


( 10 ) 

( 11 ) 


with Xl = In the limit g° = g° = gl = 0 

of no interactions within the helical modes the scaling 
dimension A 12 = < 1 is the same as for scattering 

by magnetic impuritied®. 


B. Random Disorder 

Now we switch to the case of multiple impurities on 
the edge. This perturbation is described by 

= J - V'LV'Lj)- (12 ) 

Here £,ij(x) is the (random) scattering amplitude, and 
S^ii = 0 due to the TR symmetry. We model the scatter¬ 
ing to be local along the edge and uncorrelated for the 
different pairs of helical states 

{ii 3 {x)^ki{x’)) = Wij5ik5ji6{x - x'). (13) 



FIG. 2. (color online) Phase diagram for (A®, A 2 ). Red region 
corresponds to the conducting phase for a single impurity. 
Multiple impurities are irrelevant in the blue region. The 
gray region is forbidden by positivity of matrix M. 


action the phase diagram is show in Fig. ([^ as function 
of interaction parameters. The symmetry between intra 
and inter mode forward scattering (A 4 —)■ 1 ) enhances the 
conducting phase. 
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C. Effect of two particle processes 


For random disorder, the two particle operators are 


For weak interactions, two particle processes are less 
relevant than single particle events. For sufficiently 
strong interactions, they start to compete. We analyze 
here the following processes involving two particle events 


Oii^c = ic 

0/7.s = ts 


J (15) 


,2; 

(15) 

>,2, 

(16) 


which correspond to the transfer of 2e charge and two 
particle backscattering respectively. These processes 
renormalize according to 

^ = {l-Aa)ta, (17) 


T^ii,c = J ( 20 ) 

'Dii,s = / (ix^s(a;)V'H_iV'lj,2V'L,iV'L,2, (21) 


where the ^{x)a are uncorrelated random variables with 
{^a{x)^b{x')) = Wa5abd{X — x'). These processes renor¬ 
malize acording to the RG equations 


dWa 

dl 



( 22 ) 


with Aa given by (181. For 3/4 < Fi < 1 andFi-|-F 2 > 3, 
the conducting phase remains. If Fi < 3/4 (still with 
T’l + 7^2 > 3), the process Vn^s becomes relevant under 
RG and conductance becomes non-monotonous at high 
temperatures (similar to the case of single impurity). 


with a = (s, c). Here Aq is the scaling dimension of the 

operator They are D. > 1 Helical Edge Modes 


Ac = 2 F 2 and Ag = 2Fi. (18) 

In the case of 1/2 < Fi < 1, the system is in the 
conducting phase for Fi -|- F 2 > 2. The correction to the 
conductance G = Ae^jh scales with the temperature as 


We now proceed with a more general case of v helical 
states. To calculate A, we consider the operators of the 
form 4'^ = e™'* where each vector m = (m/j,mi) 
corresponds to a different physical process. For example 


5G 


-cifJ .12 ( 


2TTaT\F+F2-2 


^ ,,2 „ (2TraT\^F-‘2 

C2fJ'12^F [ u ) 


if F 2 < 3Fi, 
if F 2 > 3Fi. 


(19) 


where Ci are non-universal parameters. If Fi < 1/2, the 
second order process Ou^s becomes relevant. The con¬ 
ductance then becomes non monotonous at large temper¬ 
atures (see Figj^ 



FIG. 3. (color online). Sketch of the conductance as function 
of temperature. For Fi < 1/2, the two particle process G/j.s 
becomes relevant, leading to a non monotonous behaviour of 
the conductance at high temperatures. In the figure we take 
F\ = K and F 2 = 1/F, which correspond to the simplest case 
of interactions just within each helical mode, discussed in Eq. 


(^^i?)z — — 6i^k+l: (23) 

describes an operator that backscatter a right mover 
in the mode fc to a left mover in the mode k + 1. Us¬ 
ing the (quadratic) action (|^, one computes the scaling 
dimension of 47^ 


(24) 

where A = A4~^A4ilCA4^\A4~F Here the absolute 
value of a matrix in the right hand side is defined as 
the absolute value of its eigenvalues. In other words, if 
H is a diagonalizable matrix A = UDU~^, then 


|H| = U\D\U-^ = U 



u-\ 


(25) 


Note that, for translationally invariant interaction we 
consider, the interaction matrices M are of Toeplitz type, 
i.e. (mf^y'r,w) 2 j = (mf^/i 3 ^)|j_j|. 

Now on we focus on the limit where the number of 
modes is large {v 1). Therefore, one can impose the 
periodic boundary condition in the mode space, without 
changing the results. In this case the interaction matri¬ 
ces are circulant, and can be easily 

diagonalizecP^. 
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We adopt a g-ology type notations and model the in¬ 
teraction by g 2 and (74 components 


(Mfw)|i_j| = vpSij + 54 ( 1 * - j|), (26) 

(Afbw)|i-j| =~52 (|*~j|)- (27) 


Here the distance dependent 54 (i) accounts for the for¬ 
ward interactions between electron densities of the same 
chirality at distance i, while g 2 {i) parameterizes the back¬ 
ward interactions of densities of opposite chiralities. The 
scaling dimension Ag, defined in Eq.(24) with m given 
by Eq.(23) is 


Ai 


1 

u 


E 




1 — G{k) cos{2'nMlv) 
\/l - G2(fc) 


(28) 


The function G{k), is determined by the interaction pa¬ 
rameters ( 72,4 


VF + 54 («) 


(29) 


Here g2,4(fc) = cos(27rjfc/:/)(72,4(j) is the cosine 

transform of ( 72 , 4 (?')- The condition of |G'(A:)| < 1 fol¬ 
lows from positivity of matrix Ai. With the scaling of 
disorder operators at hand we can analyze their behavior 
under renormalization. 

We now focus on a finite range interaction. One can 
easily show that the scattering processes between distant 
modes are less effective for the localization than backscat- 
tering between close ones. Eor the model of isotropic in¬ 
teraction g 4 {r) = g 2 {r) = g exp{—r^/R"^) the scattering 
between the distant modes is irrelevant for 
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1 R 
2vf d 


(30) 


Here we assumed that \i — j\= 1 and g/vp ^ 1- 

The scattering between the nearest modes (£ = 1) im¬ 
poses more stringent conditions on the interaction con¬ 
stants A^, as shown in Figj^ In particular, the conduct¬ 
ing phase is stable only for the nearly symmetric inter¬ 
action. For the fixed values of interaction strength the 
localization is enhanced by increasing the interaction ra¬ 
dius. In other words, strong and short range interaction 
most efficiently drives the system towards the conducting 
phase. 


III. SUMMARY 


To summarize, we have studied the localization of the 
edge modes in TIs with TR symmetry. We find that 
a combination of TR symmetry and zero bias anomaly 
changes the scaling dimensions of scattering operators. 
This notably affects the phase diagram. For a sufficiently 
strong values of interaction the zero temperature fixed 



40 


2 gs^ 4 


FIG. 4. (color online). (( 72 , 54 ) phase diagram 3 > 1 , fi¬ 
nite ranges interaction 52 , 4 (a) = 52,4 exp (—(^)^). Panels 
correspond to the different values of the interaction radius R. 
Color code is the same as in Fig[^ 


are stable against TR disorder. This holds also for the 
even fillings, where the non-interacting system is equiv¬ 
alent to a trivial insulator. We have analyzed the prob¬ 
lem in several limiting cases, for the single impurity and 
random disorder, short and long range interaction, for 
a variety of filling fractions f. We have computed the 
boundaries of the conducting phase in all these cases. 
For intermediate values of interaction electric conductiv¬ 
ity is a non-monotonous function of temperature, due to 
interplay of single and two electron scattering processes. 
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